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We present experimentally derived potential curves and spin-orbit interaction lunctions for the 
strongly perturbed A^E^ and fe^llu states ol the cesium dimer. The results are based on data Irom 
several sources. Laser-induced fluorescence Fourier transform spectroscopy (LIF FTS) was used some 
time ago in the Laboratoire Aime Cotton primarily to study the X^T.^ state. More recent work 
at Tsinghua University provides information from moderate resolution spectroscopy on the lowest 
levels of the b^Ti^^ states as well as additional high resolution data. From Innsbruck University, 
we have precision data obtained with cold Cs2 molecules. Recent data from Temple University was 
obtained using the optical-optical double resonance polarization spectroscopy technique, and finally, 
a group at the University of Latvia has added additional LIF FTS data. In the Hamiltonian matrix, 
we have used analytic potentials (the Expanded Morse Oscillator form) with both finite-difference 
(FD) coupled-channels and discrete variable representation (DVR) calculations of the term values. 
Fitted diagonal and off-diagonal spin-orbit functions are obtained and compared with ab initio 
results from Temple and Moscow State universities. 
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The lowest electronically excited states of alkali dimers 
for long have been of interest as gateway or "window" 
states for the excitation of higher singlet or triplet lev- 
els [l|. Quite recently, various low excited states of al- 
kali dimers have also been used as intermediaries in the 
production of ultracold molecules, as for RbCs LiCs 
ia], KRb ti-tS], NaCs 3> d and Csa ^-JA]- Al- 
though molecules with electric dipole moments have at- 
tracted the most interest, homonuclear species, such as 
Cs2, also offer the interesting conceptual challenge of 
non-spherically symmetric particles in a condensate . 
Furthermore the cold atoms can be prepared with fewer 



2 



lasers than required for dual species systems. Cold Cs2 
molecules may also be useful in experiments designed to 
be sensitive to the electron:proton mass ratio, as dis- 
cussed in [r^. 

Various new possibilities inherent in cold molecules 
have prompted efforts to obtain and analyze the spectra 
of alkali dimer molecules. The lowest excited states above 
the +^ S hmit, namely the A-^E^^ and b^H^u) states 
(where the ungerade designation applies only to homonu- 
clear species), would seem to be especially interesting in 
this regard. However, transitions from the ground state 
may require frequencies for which the availability of laser 
sources is limited. Furthermore, spin-orbit perturbative 
interactions between the A and b state levels complicate 
the spectra. As a route to cold ground state molecules 
from weakly bound molecules that are formed on a Fes- 
hbach resonance (FR), the A and b states of homonu- 
clear species do not have the special advantage of possible 
triplet to singlet transfer which is possible in heteronu- 
clear species 17| which have no gerade/ ungerade selection 
rules. In view of these difficulties, it is noteworthy that, 
to date, only for Cs2 have the A and b states been used 
as intermediates in the production of cold molecules [TTI - 
The spin-orbit perturbations are somewhat more 
severe in species containing Rb or Cs than in the A and 
b states of other alkali dimers. Nevertheless, with Cs2, 
Refs. [TTI - [l3 | used the X^^'^ component of a weakly 
bound molecule formed on a FR to excite to the A^Ti^ 
component of an upper state level. Two two-photon steps 
were used to efficiently connect the FR molecular state, 
via X{v = 73,J ^ 2) to X{v = 0, J = 0). 

Recently, the methods available for analyzing and 
modeling spectroscopic data on highly perturbed states 
have been extended [l8l - [27j As is evident from the dis- 
cussion below, each effort to perform a "global analy- 
sis" of all available data on these two electronic states 
raises additional questions about what Hamiltonian ele- 
ments and functions are required, and what data is re- 
quired to determine the parameters so as to achieve a 
fit with residuals comparable to experimental uncertain- 
ties. There has been impresssive progress in data acqui- 
sition and analysis for the heteronuclear A and b states 
from work on NaK [28l- [33 to more recent work on NaRb 
[Hill, NaCs [13, KCslllIii, and RbCs For 
NaRb, NaCs, and KCs, vibrational assignments of both 
states have been reliably determined, and the perturba- 
tive interactions have been modeled to an accuracy of 
0.01 cm^^ or better. Higher-order spin-orbit effects were 
included in [23| and in t25i|, and there was quite good 
agreement between empirically extracted and ab initio 
potentials. Curiously, for one alkali dimer recently used 
for cold molecule work, namely KRb 0, Q , spectroscopic 
data on low levels of the A^E+ and 6^H states appear to 
be lacking, and optical state transfer processes leading to 
cold KRb molecules have used other pathways. 

There have also been a succession of studies of the 
A and b states of the lighter homonuclear alkali dimer 
species Lis Naj [H, ^2 0, M \MM 



and Rb2 |47|, |48|. For Na2 |22|, data now extend al- 
most continuously from the lowest vibrational levels to 
the atomic limit. 

One question that arises when considering this sum- 
mary of spectroscopic data on the A and b states of Cs2 
is why hyperfine structure (hfs) is not reported. Earlier 
work [37| on the A and b states of Na2, using a molec- 
ular beam, and accessing levels with significant '^H2u 
character, showed hfs structure spanning several hun- 
dred MHz. The atomic hfs interaction of ^^"^Cs is larger 
than that of ^^Na. Two of the techniques in the present 
work, namely optical-optical double resonance polariza- 
tion spectroscopy, and spectroscopy with cold ground 
state Cs2 molecules, have very small inherent linewidths. 
Nevertheless, the transitions observed here are not ex- 
pected to exhibit hyperfine structure, as we discuss in 
Section VI below. 

This paper presents a detailed model of the energy 
level structure of the A and b states of Cs2. In section II, 
we describe experimental techniques used to obtain high 
resolution spectra. Section HI describes the model used 
to fit the data. Section IV discusses presents ab initio 
results for potentials and spin-orbit functions, calculated 
from two different approaches. Section V discusses ra- 
diative properties of the A ^ b complex, including cer- 
tain anomalies occurring for mixed singlet-triplet upper 
states. Section VI discusses the possibilities for observ- 
ing hyperfine structure, while section VII gives a short 
summary and conclusion. 



II. THE EXPERIMENTAL DATA 

This report combines experimental data from several 
sources. A substantial part comes from Fourier trans- 
form spectroscopy (FTS) performed in the Laboratoire 
Aimee Cotton (LAC), Orsay, France, in connection with 
the study of the Cs2 X^T,'^ ground state as presented in 
(49I [soj . More recently, FTS data has been obtained at 
the University of Latvia in Riga. Additional high resolu- 
tion data were obtained in Tsinghua University, Beijing, 
using fluorescence excitation spectroscopy. Another valu- 
able component of the data consists of fluorescence lines 
from the 2^Aig state observed at Tsinghua Uni- 

versity, which provided information on the lowest vibra- 
tional levels of the 6^Ho„ state, of both parities, ±(— l)"'. 
An ultra-high resolution data set, limited in absolute ac- 
curacy by wavelength meter calibration, was obtained 
by excitation of cold Cs2 molecules at the University of 
Innsbruck, covering a limited range of energies. Finally, 
because of a substantial gap in the above data on higher 
A/b levels, the Lyyra and Huennekens groups, working 
at Temple University, employed optical-optical double 
resonance (OODR) polarization spectroscopy, with col- 
lisional orientation transfer as in [48[, to partially fill in 
this gap. Unlike the other data sets, the work at Tem- 
ple has not been published previously and therefore the 
discussion is more detailed. 
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A. Data from Fourier transform spectroscopy 

Many of the term values used in the present work are 
the upper state levels from transitions observed in LAC 
and reported in [i^ and [s^. Cs atoms in a heat pipe 
oven at 600K were excited, in the 1987 study by an argon- 
ion laser, and in the 2002 work by a titanium-doped sap- 
phire (Ti:Sa) laser, pumped by an argon-ion laser. The 
fluorescence in the backward direction was collected and 
sent into a Fourier spectrometer with a 2m optical path 
length difference. Wave numbers were calibrated by ref- 
erence to a xenon atomic transition near 3.5 /im, with 
absolute uncertainties between 0.001 and 0.003 cm"^. 
For [50], 16,900 transitions were measured from 113 main 
fluorescence series and 348 series including rotational re- 
laxation. The data are no longer available on the LAC 
website, but are posted in the EPAPS file associated with 
this article fs^ ]. 

At the University of Latvia, FTS data on the A and 
b states of Cs2 have been obtained recently in the pro- 
cess of analyzing the A/b — >■ X laser-induced fluorescence 
(LIF) spectra of KCs and RbCs ^] molecules. In 

these spectra, along with transitions in KCs and RbCs, 
the A/b — >■ X LIF spectra of Cs2 molecules were present 
as well. In the experiments, the molecules were produced 
in a linear heat-pipe at about 565K temperature and ex- 
cited by diode lasers centered at 980nm, 1020nm and 
1050nm. The backward LIF spectra have been recorded 
by Fourier transform spectrometer IFS-125HR with 0.03 
cm~^ resolution. The transition frequencies obtained, 
with ca. 0.003 - 0.004 cm"^ accuracy, have been added 
to the X state term values of Cs2 calculated with the 
X state potential of [s^l, corrected by jS^j. The esti- 
mated overall uncertainty of term values, accounting for 
Doppler broadening, is about 0.01 cm~^. 

The relative LIF intensity distributions have been de- 
termined in sufficiently lengthy progressions, accounting 
for the spectral sensitivity of the InGaAs detector as de- 
scribed in Ref. [l^. For comparison with theory, the 
intensities were averaged over P, R doublet components 
and corrected for the spectral sensitivity of the detector. 

B. Fluorescence excitation spectroscopy 

This approach, employed at Tsinghua University, in- 
volved two-photon excitation to a higher state {^lig or 
^Ag) from levels of the X^Yi'^ state via levels of the 
A^E^ state, which was monitored by observing decay 
(direct, or collision-induced) to the a'^E^ state. After 
the double resonance signal was found and optimized, 
the same S^Hgiv, J ^ J' + 1) [or 2^Ag{v, J ^ J' + 1)] 
or 3^Ilg{v, J = J' - 1) [or 2^Ag{v, J ^ J' - 1)] upper 
level was excited from A^Y.+ {v', J' + 2) or A^T.+ [v', J' - 
2) intermediate levels. Thus the term values of the 
A^ll+{v', J'+2) or AiE+(i>', J'-2) levels can be obtained 
by adding the term values of the X'^Y.+ {v", J" = J' + 1, 
or J' + 3) or X^E+{v", J" = J' - 1, or J' - 3) ground 



levels [50| to the pump laser frequency with an accuracy 
of -0.003 cm-i. 



C. Spectroscopic observations on cold Cs2 
molecules 

As part of the effort at the University of Innsbruck 
to produce ultracold Cs2 molecules in the rovibronic 
ground state near quantum degeneracy (ill , [l^ , weakly 
bound molecules were formed by association of atoms on 
a Feshbach resonance in a high phase-space density sam- 
ple. Subsequently, two STIRAP transitions were used 
to transfer the population to X{v = 0). A first pair of 
lasers linked the Feshbach level to the intermediate level, 
X{v = 73, J = or J = 2), via an A/b excited state 
level near 12,500 cm^^. Unfortunately, in the present 
study there were insufficient data points near the levels 
around 12,500 cm~^, so these term values could not be 
accurately fit in our analysis. 

The second two-photon transition starting from the 
X{v — 73) level involved an A/b intermediate level near 
10,000 cm~^. Spectrosco py on the first transition, as 
described in detail in Ref. [13| , was carried out by irradi- 
ating the Feshbach molecules for a given time with light 
near 1126 nm, stepping the laser frequency with each 
cycle of the experiment. For spectroscopy near 10,000 
cm~^, molecules were first coherently transferred by STI- 
RAP to X{v — 73, J — 2), and then irradiated with 
light near 1350 nm, as described in Ref. [IJ]. For de- 
tection, the molecules were transferred back to the orig- 
inal weakly bound state, dissociated, and the resulting 
atoms were imaged by absorption imaging. For STI- 
RAP, the lasers were stabilized to reach both a short-term 
linewidth and long-term stability on the kHz level by 
locking to optical resonators and additional referencing to 
a stabilized optical frequency comb. However, waveme- 
ter calibration was estimated to be 0.01 cm~^, presenting 
a limitation on the absolute frequency accuracy of these 
transitions. Ultimately, the radiative linewidth of per- 
haps 1 MHz would limit the resolution. 

Approximately lA A/b levels excited from X state lev- 
els near w=73 are included in the data set analyzed in the 
present work. These term values are listed in Table 1 of 
|12| in connection with an earlier version of the present 
analysis. 



D. Fluorescence from higher- lying states 

These data were obtained at Tsinghua University, and 
are taken from Ref. [5l|. Two successive diode lasers, 
each of 5 MHz linewidth, were used to excite A^Yj'^ ^ 
X^E+ and 2^Aig -f- A^E+ transitions. Fluorescence 
from 2^Aig — > &'^Ho„ was dispersed with a 0.85 m double 
grating Spex 1404 monochromator. The accuracy of the 
fluorescence lines was estimated to be about 1.5 cm~^, 
with some better and others worse. 
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E. Data from OODR polarization spectoscopy, 
facilitated by collisional molecular orientation 
transfer 

V-type optical-optical double-resonance (OODR) po- 
larization spectroscopy was performed at Temple Uni- 
versity to obtain data on mixed A/b levels in a region of 
energy not studied previously. Among different types of 
polarization spectroscopy |54l - [57| , this technique is well 
suited to observe selected transitions in the congested 
spectra of Cs2 since in the single-laser based sub-Doppler 
experiments, the assignment of spectral lines is not triv- 
ial. Furthermore, collisional transfer of orientation pro- 
duces a series of well-separated spectral lines associated 
with nearby levels that have been oriented and polarized. 
Another possible approach might be optical triple reso- 
nance, as used for Lis [111, Na2 [13, and K2 [ij [ist- 
This approach would overcome the unfavorable Franck- 
Condon factors mentioned below, but would not yield the 
range of J values obtained in the present work, since opti- 
mization of signal to noise through cooling of the sample 
molecules would probably be necessary. 

As in previous implementations of this technique (see 
(isj and [l^l and references therein), a circularly polar- 
ized pump laser beam, tuned to a chosen transition, cre- 
ates an anisotropic distribution of magnetic sublevel pop- 
ulations (net orientation) in one or more ground state 
rovibrational levels. The counterpropagating linearly po- 
larized probe laser beam can be considered to be made up 
of equal parts left and right circular polarization. When 
the probe frequency is tuned to a transition sharing ei- 
ther level with the pump transition, the orientation cre- 
ated by the pump causes the two circularly polarized 
probe components to experience different absorption co- 
efficients (circular dichroism) and different indices of re- 
fraction (circular birefringence) . Consequently, when the 
probe beam exits the oven, the two components no longer 
sum to the initial linear polarization, but rather the beam 
has a slight elliptical polarization and some fraction of it 
is transmitted through a final crossed linear polarizer be- 
fore reaching the detector. Figure [T] shows schematically 
the transitions in this technique. More information on 
the C^Htj state is available from the potentials shown in 
Fig. [S] A representation of the apparatus is shown in 
[22| . The experimental setup used here was similar to 
the one used previously for Rbs fis']. 

Collisional lines help to expand the data field but can 
only be observed if orientation is transferred from the 
ground state level labeled by the pump laser to a neigh- 
boring level. A quantitative study of the transfer of pop- 
ulation and orientation in collisions of NaK molecules 
with argon and pot assium atoms has been carried out 
and reported in [SSj and [s^. These references discuss 
the transfer of orientation following the discussion of Ref. 
[60| , and also present a complete analysis of polarization 
spectroscopy line shapes, extending the textbook presen- 
tation in [6l!|. 

In our experiments, the Cs2 metal was loaded at the 



center of a five-arm heat pipe oven, which was kept 
at a temperature of 550 K, with 3 Torr Argon buffer 
gas. A CR-699-29 tunable single mode laser with Ki- 
ton Red 620 dye was used as the circularly polarized 
pump laser, tuned to known [6^] Cs2 C^n„(w', J' ± 1) ^ 
Xii;+(u", J") transitions. A CR-899-29 Ti-Sapphire sin- 
gle mode laser served as the linearly polarized probe laser. 
The two lasers were counter-propagating and overlapped 
in the center of the heatpipe with a crossing angle that 
was minimized so as to maximize the overlap area. The 
power of the pump (probe) laser was lOOmW (20mW) in 
front of the heatpipe. The spot sizes (defined as a ra- 
dius at 1/e^ intensity) of the pump (1.2 mm) and probe 
(0.9mm) laser beams in the interaction region were mea- 
sured with a razor blade technique [6^ . Two lenses were 
placed in the path of each laser beam for collimation in- 
side the heatpipe. 

A Glan Thompson linear polarizer and a A/4 Babinet- 
Soleil compensator from Karl Lambrecht Corp. [63j were 
used to make the pump beam circularly polarized. Two 
linear polarizers in the probe laser beam path (one before 
and one after the heat pipe oven) had relative transmis- 
sion axes at a 90 degree angle for maximum extinction of 
the probe laser when it was off resonance. A power meter 
(Coherent LabMax TO) with a sensitive detector (LM-2) 
placed after the second polarizer was used to make fine 
adjustments of the orientation of the second polarizer to 
create an optimal extinction ratio of 10~^ (and then re- 
moved for data acquisition). The observed ratio includes 
effects of window birefringence and circular dichroism. 
A photomultiplier tube (Hamamatsu R636-10) detected 
the polarization signal. Its output was amplified using 
a lock-in amplifier (Stanford Research Systems SR 850) 
for phase sensitive detection. The pump laser was modu- 
lated at a frequency of 980 Hz using a mechanical chopper 
(Stanford Research Systems SR 540). 

The pump laser was calibrated using an iodine atlas 
[gsI [66| . with estimated accuracy of ^0.004 cm~^. The 
Ti:Sapphire laser was calibrated by comparing optogal- 
vanic spectroscopy signals from a uranium lamp to line 
positions in the uranium atlas (67j when it scanned above 
11000 cm^^. The listed wavenumbers in the Uranium 
atlas are accurate to ±0.003 cm~^. A Burleigh 1600 
wavemeter and BOMEM DA8 FTIR were used to cal- 
ibrate the Ti:Sapphire frequency when it scanned below 
11000 cm~^. The calibrated wavenumbers of the exper- 
imental peaks were added to the term values of their 
corresponding ground state levels, calculated from the 
Dunham parameters of [sol], to get the A/b state term 
values. 

Because the Franck-Condon factors for transitions of 
interest were weaker by a factor of 100 in Cs2 than for 
the Rb2 lines in [48] , the signal to noise ratio was smaller 
and fewer collisional satellites were observed than with 
Rb2 (AJ < 12 rather than < 58). Also, due to spectral 
congestion in Cs2, the pump transitions were difficult to 
isolate. Some of the pump transitions were confirmed by 
analyzing the resolved fluorescence recorded by a SPEX 
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model 1404 0.85 m dual grating monochromator. Other 
pump transitions were confirmed by scanning the pump 
laser while the probe laser was tuned to a previously as- 
signed coUisional peak. 

26 different C X transitions (ground state 3 < w < 8 
and 52 < J < 121) were confirmed and 20 of them served 
as pump transitions in the experiment presented in this 
work. These transitions are listed in the EPAPS file [l^l • 

Figure [2] shows one OODR polarization spectroscopy 
data scan. For homonuclear molecules, the selection rule 
for collisions is A J = even due to the symmetry proper- 
ties of the nuclear spin wavefunctions (for states that are 
not parity doubled). However, in the case shown in Fig. 
2, the C^IIk and -^^5]+ rotational spacings are such that 
two pump transitions Cin„(7,73) - XiE+(6, 74) and 
C^Iiu{1, 72) - XiE+(6, 73) overlap within their Doppler 
widths and are therefore simultaneously excited. Thus 
Fig. 2 shows two separate collisional series (both charac- 
terized by A J = even) starting from labeled ground state 
levels (6,74) and (6,73). In this scan, the off-resonance 
baseline is artificially created by a non-zero offset added 
by the lock-in amplifier. The P branch lines are positive- 
going relative to this offset, and the R branch lines are 
negative-going here. The lock-in was used here in the 
"R" mode, in which the reported signal is the square 
root of the sum of the squares of the "in phase" and "in 
quadrature" components, and hence must always be pos- 
itive. Thus if the R branch lines are particularly strong, 
they will reflect off the zero signal baseline, appearing as 
a central peak with minima on either side, as shown in 
Figs. 3b and 3c. The fact that R branch and P branch 
lines have opposite signs (as always occurs if an offset 
is used or if only the "in phase" channel is recorded) is 
helpful for making assignments. In [48], the lock- in "R" 
mode was also used, but without an offset, so P and 
R lines both appeared positive going in the polarization 
scans reported in that work. 



F. X state term values 




A^x; (v',j"+i) 

A^X;(v',J"-1) 



FIG. 1. A schematic of the transitions used for OODR polar- 
ization spectroscopy. 
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Recently Coxon and Hajigeorgiou 53| have reanalyzed 
the data of Amiot and Dulieu [50] on Cs2 A — X transi- 
tions to make a direct fit to a potential of the X'^Yj'^ state, 
and hence to obtain improved X state term values. It has 
come to be recognized, as asserted in [gl], that using a 
direct fit to an analytic potential is more precise than the 
semiclassical Dunham procedure. The differences in this 
case are less than 0.01 cm~^, but nevertheless significant 
compared to residuals from least squares fits to A—h state 
term values presented below. Supplementary data asso- 
ciated with [5 31 give X state term values for transitions 
in the data set of (HO]. For other transitions used in this 
study, we have calculated X state term values from the 
potential of [5^. In doing this, we have used the discrete 
variable representation (DVR) described below, to obtain 
eigenvalues at high v and high J levels without the use 
of high-order Dunham coefficients. When compared with 



FIG. 2. OODR polarization spectroscopy signals. For 
this scan, the pump laser was set to 15898.5535 cm~^, si- 
multaneously labelling the ground states of the transitions 
C'n„(7, 73) - X'-T,+ (6,74) and C^Tlu{7, 72) - X^Y.+ (6,73). 
In this scan, P branch lines are always positive going, and R 
branch lines are negative-going, except for the strong lines, 
ii(73) and _R(74), that reflect off the baseline as described 
in the text. Note that transitions from the J=74 pumped 
level are relatively strongest. Those associated with other 
lower state rotational quantum numbers are due to collisional 
transfer of orientation. 



X state term values listed by [53| , these values agreed to 
within 4 x 10"'' cm~^. 
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FIG. 3. OODR polarization spectroscopy line shapes for the 
particular conditions used in the present experiments. In each 
panel, there is a positive-going P line and an R line. In (a), 
the R line is weak and entirely negative- going. In (b) and (c), 
the R line is stronger and reflects off the zero signal baseline, 
producing a double minimum line shape. The vertical scale 
varies between parts (a)-(c), as can be seen by comparison 
with the previous figure. 

G. Summary of the data 

Figure|3]sunimarizes the experimental data used in this 
study. It will be noted that the data from Temple Uni- 
versity partially fill in a gap, and that (high-resolution) 
data from Tsinghua and Riga augment the data available 
near the minimum of the A state, as shown more clearly 
in Fig. [To] The various gaps in the data are associated 
with regions of low Franck-Condon factors (see Sec. |V]) 
and unavailability of lasers of suitable wavelength. Data 
from Innsbruck, at low J values, extend up to 12,600 
cm^^ above the minimum of the X state. However, be- 
cause data in this region was very limited, we were not 
able to extend the least squares fits to the data beyond 
a reduced energy Ered = E - 0.0091J( J -f 1) of 12,400 
cm~^. 

When working with data from different sources, there 
is always a question of relative calibration errors. The 
average residuals (shown in Table U below) indicate that 
in this case differential calibration errors are smaller than 
the rms residuals from the fitting process. There were 
few cases of overlapping data to pose a more severe test 
for calibration differentials. Six term values measured 
by LAC and Tsinghua displayed an average difference of 
0.0025 cm-i. 



III. ANALYSIS OF THE DATA 

Because the data are relatively sparse in certain re- 
gions, the analysis has required special procedures and 




FIG. 4. Term value data used for this study, (a) Higher 
resolution FTS, OODR and cold molecule spectroscopy term 
values, (b) and (c) display lower resolution (monochromator) 
data from Tsinghua University, on e and / parity levels, re- 
spectively. Energies in this plot are relative to the minimum 
of the X^E^ ground state. The Innsbruck data above 12,500 
cm^^ were actually not used in the fit. However, all other data 
shown here fall within the limit Ered = E- 0.0091J( J -I- 1) < 
12400 cm~^, and are included in the least squares fit. Please 
see Fig. [10] for a clearer picture of the extent of data near the 
minimum of the A state, for example. 

furthermore the quality of the fit is not quite comparable 
to that of Ref. [2^, for example. Many of the methods 
are similar to those used in previous studies [1^ HI] , and 
therefore will be summarized briefly, with attention to 
the special procedures mandated by the available data 
set. It is important to note also that for this relatively 
heavy alkali dimer, the second-order spin-orbit (S02) ef- 
fects become significant, and will require special atten- 
tion below. 

Hamiltonian elements. The molecular Hamiltonian 
M 

H = Hbo + Hk + Hso + Hrot, (1) 

includes the Born-Oppenheimer potentials Hbo, radial 
kinetic energy Hk, nuclear rotation Hrot, and spin-orbit 
interaction Hso- Hyperfine structure has not yet been 
observed, but in Sec. I VII we attempt to estimate the 
magnitude of the hyperfine structure in the relevant 
states of Cs2 . 

The observed term values consist of levels of b^^"^^ 
below v=Q of the A^Ti^ state, together with many levels 
that can be identified as mixtures of ^^S^ levels and 
6'^nJ„ levels. In the first category are levels of both parity 
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(— 1)'^ (e symmetry) and parity —(—1)'^ (/ symmetry). 
Overall, only a few levels were found to have significant 
■^IIi character and no levels with significant ^112 character 
were identified. In view of the very limited information 
on the "^IIi levels, our '^IIo potential contains the spin- 
orbit function implicitly. When '^IIi and ^112 levels are 
considered, the matrix elements of Hbo + Hso + Hrot + 
Hso2 are 



{'K\H\'K) 




X 4 


- 2)B + C/l°^ 


(3n+|ij|3n+) 


- Ko + 


{x^ 


- 2)B + 




= + 


{x^ 


-2)B- 


1 TTS02 

r <^bOf 


fni„|ir|3ni„> 


= u+ + 


{x^ 


-2)B- 




(3n2„|F|3n2„) 


= K + 


[x - 


- 2)B- 


H'bio+C 


CEj|iJ|3n+) 


= -V2ao 







(3n 



ou\H 
iu\H 



/iv+ 



~,b21 



(2) 



2m) 



= -V2xB 



3n 



= -^J2{x-2)B 
= -V2xB'q 



where x — J( J+1), and — H, where is the trans- 
pose of H. In the above, the potentials, Ua, and U^, as 
well as the spin-orbit functions, ^^"^q (off-diagonal) and 
CbiO' '^621 (on-diagonal); and also B = ?i^/2/ir^ are func- 
tions of internuclear distance r {fj, is the reduced mass). 
^610 = 1 — 0"'" interval of the b state, and ^^21 

is the n — 2 — 1 interval. Both and ^^21 functions 
include the S02 effect implicitly. The parameter 77 arises 
from second-order spin-orbit and electronic-rotational ef- 
fects, as discussed in (2^ . 

The second-order spin-orbit (S02) terms in the above 
Hamiltonian elements have been estimated by ab ini- 
tio calculations using the present quasi-relativistic ECP- 
CPP-CI procedures described in Sec. IIV Al These func- 
tions are plotted in Fig. [5}d and given numerically in the 
EPAPS [Hi file. The S02 functions plotted in Fig. Eb 
are: 



Ua - Uj 
214^0-/ 



bO 



j e (2)3n„ 



J e (1 - 3)3I]^ 



bl 



(3) 
(4) 
(5) 
(6) 



J e (i-3)3i]J;(i-2)in„ 

From Fig. [5)d, it is evident that C/go/ large and positive 
for r < 5.2 A, and becomes negative for large r. This 
behavior can be understood from the b^H^^ — a'^Sj spin- 
orbit coupling function, which shifts the b^IlQ^ energies 



Cs(6^S+5^D) 




R(A) 



FIG. 5. (a) Potentials for ungerade states of Cs2 dissoci- 
ating to the three lowest atomic limits, as calculated from 
ECP-CPP-CI methods, (b) Second-order spin-orbit shifts 



for A^T,i 



and b IIiuc calculated from these potential 



functions and the spin-orbit coupling functions, using Eqs. [S] 
llElandEl 



potential lies just 

below the potential in this region. At larger values 



upward at small r, because the a" 



of r, the b'^Hg^ potential approaches various ungerade 
potentials from below, as seen in Fig. [5^, and thus the 
S02 shift is negative. 

Analytic potential and spin-orbit forms. Al- 
though an ideal deperturbation analysis should yield the 
identical Uj^ and Uj^ PECs (potential energy curves) for 
the triplet state, the experimental term values assigned 
to the e(-l-) and /(— ) parity levels were involved in the 
fitting process separately. The empirical PECs of the in- 
teracting A"^E+ and ^'^Hg^ diabatic states Ua, Uj^ were 



represented analytical^ 
tor (EMO) function [tC 



by the Expanded Morse Oscilla- 



(7) 
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converging to the appropriate dissociation limit, namely: 
Tf = Td^s - Df for the the singlet ^-state and = 
Tms - £.cs - Df for the = sub-state of the triplet b- 
state. Here ^g'^ = [-Bg^Pa/^ "-^e^Pi/J/S = 184.6797 cm"! 
is the empirical spin-orbit constant of the Cs atom in the 
&P state [lil while Tdis = D^+Eq2p-Eq2s is the energy 
of the center of gravity of the Cs(6^P) doublet (without 
hfs splitting) relative to the minimum of the X state po- 
tential. From |72i Eqip-Eq2s = 11547.6275 cm^^. The 
dissociation energy of the ground state (for the hfs center- 
of-gravity) = 3650.0299(14) cm"! from Ref.[ll|, tak- 
ing into account a correction of -0.0022 cm~^ for the en- 
ergies oiX{v = n- 73) levels, following [H]. Thus Td^s 
= 15197.6574 cm'^. Note that because of the [71°^ and 
fJ^Q ^ terms, the effective potential minima differ from the 
fitted Te values. 

Both off-diagonal ^^°^q and on-diagonal ^^^qj Cmi 
pirical spin-orbit coupling functions were approximated 
by the Hulburt-Hirschfelder (HH) potential [73| 

r M = Cc°s - DT [2e-- - e-2-[l + cx\l + bx)]] ;(8) 
X ^a{r/rf - 1) 

determined by the five fitted parameters Dl° , r|°, a, b 
and c. 

The initial EMO and HH parameters of the relevant 
PECs and SO functions were estimated using the present 
quasi-relativistic ECP-CPP-CI electronic structure cal- 
culation. The refined parameters of the PECs and SO 
coupling functions f^°{,Q and ^^^q were determined iter- 
atively using a weighted nonlinear least-square fitting 
(NLSF) procedure. As there were not data on 6'^H2„, 
the initial ah initio estimate for ^^21 used in the least 
squares fits. 

Fitting procedures. The spin-orbit coupling func- 
tion between the and &^HJ„ states of Cs2 is larger 
than the vibrational intervals, thus requiring careful 
treatment of perturbative interactions. Initially, the least 
squares fitting was performed with a 2x2 Hamiltonian 
matrix, including simply the £7 = 0^ manifolds and the 
off-diagonal spin-orbit coupling term, ^^°{,o- 2 chan- 
nel fit was crucial in determining nearly final potentials. 
The details of the robust fitting procedure used are given 
elsewhere [13, HH] • Next a 3 x 3 Hamiltonian included 
•^Hi levels and ^^°q. Final fits were made with a 4 channel 
Hamiltonian and ^^21- Although ^H2 components were 
not directly observed, they affected the fitted parame- 
ters through the terms coupling with ^Hi and thus indi- 
rectly to ^IV^ ■ These multi-channel calculations utilized 
both finite-difference (FD) and discrete variable repre- 
sentation (DVR) approaches. For additional details on 
the FD and DVR numerical methods, please see [2^ and 
(isj , respectively. It was determined that if identical po- 
tential and spin-orbit functions are used, the eigenvalues 
obtained from the FD and DVR methods are in excellent 
agreement. 

Fitted parameters, potentials, spin-orbit func- 
tions, and term values. The e parity term value data 
was fit with 12(8) potential parameters, ai in Eq. [71 for 
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FIG. 6. Residuals from least squares fits, (a) shows residuals 
from fits to the high resolution data from all four sources; (b) 
shows residuals from the fits to low resolution (monochroma- 
tor) data. 



TABLE I. Residuals from the least squares fit to Cs2 A/b state 
data. Np — number of data points. Rms. = Rms residual in 
cm~^; Ave. — Average residual (observed minus calculated 
term value), also in cm~^. Tsinghua(LR) denotes the low 
resolution (monochromator) data from Tsinghua University 
([sH): only e parity levels were included in this fit. The exper- 
imental errors are approximately 1.5 cm^^ for Tsinghua(LR) 
and 0.01 cm~^ for the other data sets. 



Source 


Np Rms. 


Ave. 


LAC [50, 53] 


340 0.0154 0.0072 


LAC [49| 


55 0.0198 


-0.0045 


Tsinghua 


58 0.0086 


0.0028 


Innsbruck 


14 0.0075 


0.0010 


Temple 


161 0.0114 


-0.0006 


Riga 


75 0.0098 


0.0041 


Tsinghua(LR)(e) 194 1.554 


-0.437 



the 6'^ng^(A^S+) state, plus re and Te for each state. 
The HH forms (Eq. © used for ^J^g and required 5 
fit parameters each, while the parameters for ^^21 were 
provided by ab initio calculations (ECPl), The range of 
the fitted data extended from 7,982 cm"^ (w=0, J=12 of 
6'^Hou-i-) to 12812 cm~^, J = 326. But a more meaningful 
measure is the reduced energy, Eres = i?— 0.0091J( J-f-l), 
which extended only to 12390 cm~^. By contrast, the 
dissociation limits of the b^Hou and states lie at 

15012.98 cm-i and 15197.66 cm-1, respectively. The 
S02 shift terms, calculated by Eqs.®, Q, © and given 
numerically in the EPAPS files [52], were added to the 
A^J:+ and b^U^^ potentials. 

Because of the relatively large uncertainty in the S02 
terms, data for b^Hg^ (/ parity) were not included in the 
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fit. 

The quality of tlie final fit to the e parity data is in- 
dicated by the plotted residuals in Fig. |6l and displayed 
for the individual data sets in Table HI Table HIl gives the 
fitted potential parameters together with results from ab 
initio ECP-CPP-CI calculations. Table Hill gives parame- 
ters for the empirical off-diagonal and diagonal spin-orbit 
functions. Figure [7] displays the empirical potentials for 
the A^T,^ and fo'^Hi^^ states over the range of the data, 
namely up to 12,400 cm~^ above the minimum of the X 
state. 

Numerical values for the potentials and spin-orbit func- 
tions, parameter listings, calculated term values and 
comparisons between calculated and observed term val- 
ues, and several plots of calculated and observed term 
values are all given in the EPAPS files (52j . 

It should be noted that in fitting the data for -f (e) par- 
ity levels from the Tsinghua monochromator (low resolu- 
tion) data, it was found that assigning certain terms to 
^n2„ rather than '^Ilg^ reduced the residuals. Ultimately, 
such assignments were not accepted because they violate 
known selection rules for pure (a) Hund's coupling case. 
The larger residuals obtained when these levels were as- 
signed to ^n^^ may be associated with the low intensities 
of some of the lines as noted in Ref. [5l|, or they may 
be associated with breakdown of Hund's case a selection 
rules due to spin-orbit effects. 

The off-diagonal spin-orbit function ^^^g connecting 
A^Ti^ and b^^^^ states is determined empirically from 
the perturbation crossings. The fitted parameters are 
given in Table lllli and the fitted function is plotted in 
FigM together with ab initio functions for which the com- 
putational methods are discussed below. The fitted ^^°{,q 
function is primarily determined by its value at the po- 
tential crossing point, denoted in this figure. It is 
satisfying that at r^,, the fitted function is indeed close 
to all three ab initio functions plotted also in this figure. 
For other values of r, there are significant differences. 

The fine structure splitting between the b'^Ho and b^Hi 
potentials could be deduced from the presently available 
data only from a few regions of avoided crossings in the 
rotational progressions. We have found four such cases, 
three of which are shown in Fig. [9] From these deviations 
in the rotational level energies, we have fit the parameters 
(see Table IIII[) that characterize the fine structure inter- 
val. The fitted ^^°q function and an ab initio result are 
plotted in Fig. |8l The minimum of this empirical fine- 
structure splitting function is approximately 139 cm~^, 
as compared with the minima of the analogous functions 
for RbCs and KCs, reported to be 71.4 cm~^ [23| and 
78.26 cm~^ [2^, respectively, both based on more com- 
plete data than is presently available for Cs2. The 
ECP-CPP-CI ab initio spin-orbit functions happen to be 
in good agreement with the empirical function. It should 
be noted that fits of comparable quality were obtained 
with an empirical ^^°g function exhibiting a minimum 
of ~ 120 cm~^, hence with a different vibrational num- 
bering for the b^Iliu manifold. The preference for the 



TABLE II. Parameters obtained from a fit to EMO poten- 
tials using Eqs.(6) and (7), together with ab initio results. In 
each case, r,.e/ ~ 5.0 A and p = 3. re values are in A. Te, 
De = Tdis — Te and uje are in cm~^, and the other parameters 
are dimensionless. i] = —0.08228330 is defined in Eq.(2). Ab 
initio ECPl-CPP-CI values for b^Ylou are obtained from cal- 
culated f/(^ni) - The (jJe values are obtained from the 
second derivative of the potentials at r = re. 



Expt. ah initio 







A 2-1^ 






Te 


4.45746 


5.32913 


4.5031 


5.3799 


Te 


7977.8479 


9587.1155 


8077.8 


9809.8 


fflo 


0.522835586 


0.4332952472 






ffll 


0.135607457 


0.0431537069 






02 


0.166158407 


-0.01827733744 






as 


-0.097477983 


0.12061765970 






04 


-0.521843942 


-0.09455975827 






05 


1.9047134070 


0.08598874346 






oe 


1.737199081 


0.42659263834 






Ot 


-7.992154379 


-0.47446387217 






OS 


0.076434068 








Og 


8.921002546 








010 


-3.390109463 








De 


7035.1427 


5610.5398 






UJe 


42.6524 


33.0055 


43.133 


31.953 



TABLE III. Parameters for the spin-orbit functions in the 
Hulburt-Hirschfelder form given in Eq. (8). The parameters 
for ^^"^,0 and i^u) ^re from the least squares fit, while those 
for ^521 s-re from ab initio calculations. In each case, = 
184.6794 cm~^. _D|° values are in cm~^, r|° are in A, and a, 6 
and c are dimensionless. 

so ^SQ ^SO 

AM) ShlO Kb1\ 

Dl° 70.87086 46.10116 37.2907 
rl° 5.968745 6.288671 6.341088 
o 2.25268 3.34271 3.13278 
b 1.221883 0.473373 0.421783 
c 0.325679 0.4050794 0.3027172 

function plotted in Fig. |8]was based on close agreement 
with ab initio results, in view of the good agreement be- 
tween empirical and ab initio SO functions in the case of 

KCs m. 

The last two columns of table HIl give ab initio ECPl- 
CPP-CI results for parameters derived from J7(A^I]+) 
and from U{b^IV^^) — ■Cbio- It is also interesting to com- 
pare results from our least squares fit to the experimental 
data with earlier calculations by Spies [t^], which were 
not published but were widely circulated among people 
working with Cs2. However, since these calculations were 
primarily relativistic, they pertain to the adiabatic po- 
tentials, as shown by the dashed lines in Fig. [T] To make 
a comparison, we have calculated adiabatic potentials 
from the experimental potentials and from the ab ini- 
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FIG. 7. Potentials for the Cs2 A^E^ and b^Ilu states obtained 
from fits to the data in this work, over the range sampled by 
the data. The b^Il2u potential (uppermost of the set of three) 
is obtained from ab initio calculations only. The dashed lines 
indicate the adiabatic potentials, obtained by diagonaliz- 
ing the diabatic A and 60 potentials with the spin-orbit cou- 
pling term, ^Abo- 



TABLE IV. Parameters fit to adiabatic 0^ potentials from 
various sources. Near r^, the lower adiabatic potential is close 
to the &^nj„ potential, while the upper adiabatic potential is 
close to the A^E^ potential. Te and We are in cm~^, while Te 
is in A. 





Expt. 


ECPl-CPP-CI 


Spies 


Te{lo'wer) 


7960.45 


8064.79 


8169.57 


re{lower) 


4.4583 


4.5049 


4.5028 


ujeilower) 


42.432 


42.904 


43.082 


Te{upper) 


9626.64 


9843.53 


9598.16 


Veiupper) 


5.2903 


5.3465 


5.3327 


ujeiupper) 


36.541 


34.066 


35.135 



tio ECPl-CPP-CI potentials (by diagonalizing the poten- 
tials plus appropriate spin-orbit functions as a function of 
r), and then we have extracted values for Te,re, and Wg, 
to compare with similar quantities fit to the potentials of 
ref. [74]. Table ITVl shows that, considering the inherent 
limitations of quantum chemistry calculations, the ear- 
lier ab initio results were in moderately good agreement 
with current experimental results and ah initio results. 

Energy level structure - e parity levels. Figure 
[To] displays term values calculated from fitted parameters 
together with the input term value data, on a reduced en- 
ergy scale to flatten the rotational structure. One point 
to note here is that the experimental data for low lev- 
els of the A state are quite sparse as compared with the 
study of the A and b states of Rb2 [4g|, for which it 
was possible to observe long rotational progressions that 
facilitated accurate fitted parameters. The top part of 



FIG. 8. Fitted and ab initio off-diagonal spin-orbit function, 
\/2 X ^AbO) coupling the A^E+ and b^^^^ states, together with 
the diagonal function, ^^jq, that gives the fine-structure split- 
ting between 6^n^„ and 6^n+„. Both MR-RAS-CI and ECP- 
CPP-CI procedures were used for the ab initio functions, as 
indicated. For the latter, the upper (lower) line represents re- 
sults from ECPl (ECP2) basis sets. The vertical line denotes 
the value r = of the potential crossing point. For ^fJoi the 
ECPl and ECP2 results agreed to within the linewidth. 
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FIG. 9. Examples of avoided crossings between f2=l and 
n = 0"*" levels, (a) shows the relatively steeper slope of $1=1 
levels in a plot of reduced energy vs. J(J -I- 1). (b) is an 
enlargement of the crossing region in (a), (c) shows another 
crossing region. In (d), a crossing between mixed $7 = 0''' 
levels lies very near to an = 1 level. Small xs denote 
calculated term values, larger filled circles are observed term 
values in these plots. 
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Fig. [To] displays a region with more ample data. Note 
that, despite the strong perturbative coupling, the effec- 
tive rotational structure clearly differentiates states that 
are primarily A^T,^, which have the least slope, from 
those that are primarily b^IlQ^, which have a slope inter- 
mediate between the former and the b^Hiu states. In this 
regard, please note in the top part of Fig. [101 that several 
of the levels observed in Innsbruck University by excita- 
tion from cold molecules in levels of the X^T,'^ state, 
appear to have primarily triplet character. This feature 
will be discussed further in Sec. [Vj in connection with 
transition amplitudes between mixed A ^ b levels and X 
state levels. 

The contrast in rotational structure (slope vs. J{J + 
1)) of A and b state levels seen in Fig. [TUj persists up to 
approximately v=AO of the A state. Higher levels, such 
as those shown in Fig. [TT] exhibit a large degree of mix- 
ing such that it becomes impossible to assign vibrational 
quantum numbers. 

Energy level structure - f parity levels: Origin 
of the e — / A-doubling effect in the b'^llon sub- 
state. As was mentioned above the main contribution 
to the ri(A)-splitting of the b^XIoM sub-state 



E, 



vJ 



E. 



vJ 



(9) 



comes from the strong spin-orbit coupling with the 
nearby sing let A'^Y.+ state, the "SO" contribution. How- 
ever, a significant contribution ("S02") comes also from 
the second-order spin-orbit shifts, which are written in 
Eq. [SI and plotted in Fig. [S)d. Theoretical values for the 
"SO" contribution are obtained by taking the difference 
between eigenvalues calculated with the fitted potentials 
with and without the off-diagonal ^^°^q coupling func- 
tion. Calculated values for the "S02" contribution are 
obtained as an expectation value from the second-order 
perturbation shift function, U§q^, shown in Fig. [SJj, and 
the rovibrational wavefunctions calculated from the Uf^Q 
potential. The results are plotted in Fig. [121 There is 
a rapid rise as the b^^IIou levels approach A{v — 0). The 
contribution of S02 terms increase more slowly between 
3.5 and 5 cm~^, as the fo'^H^^^ functions sample the S02 
shift function shown in Fig. [Sja. Because the S02 terms 
are likely to be uncertain by 5-10% (hence possibly as 
much as 0.5 cm~^), we have not attempted an empirical 
fit to the / parity data. There may also be an " intrinsic" 
e — / energy difference due to spin-spin interaction terms, 
but but we have not attempted to estimate such effects. 



IV. AB INITIO CALCULATIONS OF 
POTENTIALS AND SPIN-ORBIT FUNCTIONS 



In this section we discuss the methods used for ab initio 
calculations at Moscow State University and at Temple 
University. 
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FIG. 10. Term values near A{v)=0 (bottom), showing spar- 
sity of data, and at higher energies (top), showing somewhat 
more abundant data from various sources. Open diamonds 
indicate monochromator data, open circles denote high res- 
olution data, also from Tsinghua, closed circles denote high 
resolution data from LAC, Xs denote data from cold molecule 
spectroscopy from Innsbruck, and triangles denote data from 
Riga. The most steeply sloping eigenvalues are for f2 = 1 lev- 
els, the least steeply sloping are for levels that are primarily 
A^Tjt in character. 



A. Quasi-relativistic ECP-CPP-CI calculations (at 
Moscow State University) 

The potential energy curves, transition dipole mo- 
ments, spin-orbit and angular coupling matrix elements 
between the Cs2 electronic states converging to the low- 
est three dissociation limits were evaluated in the basis 
of the spin-averaged wavefunctions corresponding to pure 
Hund's coupling case (a). The quasi-relativistic matrix 
elements have been obtained for a wide range of inter- 
nuclear distance by using effective core pseudopotential 
(ECP) [t^. The core- valence correlation has been taken 
in account using a large scale multi-reference configu- 
ration interaction (MR-CI) method [t^ combined with 
semi-empirical core polarization potential (CPP) [77l |. 
All calculations were performed by means of the MOL- 
PRO V.2008 program package [zll. 

To include relativistic effects, the inner core shell of the 
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FIG. 11. Some of the term values (diamonds) measured 
at Temple University to partially fill a gap in earlier data. 
Smaller filled circles denoted calculated term values. Note 
that in this energy region, the 0^^ energy levels repel each 
other and do not exhibit narrow avoided crossings, as they do 
in the previous figure. 
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FIG. 12. Energy dilTerence between e and / parity states for 
levels of b^Tiou, as a function of vibrational level, v, and rota- 
tional level. Results are plotted for the 5 values of J observed 
in [12], namely J=12, 24, 44, 78 and 100. These J values are 
translated into fe=0,4, and the horizontal axis is u + fe/S, so as 
to spread out the data on different J levels. The theoretical 
S02 contribution is obtained from the second-order spin-orbit 
coupling functions and the potential energy differences shown 
in Fig. \S\ The theoretical " SO" contribution is the difference 
between the calculated energy with and without ^^bo- "Total" 
denotes the sum of the SO and the S02 terms. Dots denote 
theoretical values, x symbols denote experimental values. 



Cs atom has been replaced by spin-orbit averaged non- 
empirical small core 9-electrons ECP, leaving 18 outer- 
core and valence electrons of the cesium dimer for ex- 
plicit correlation treatment. In order to test the ECP 
basis set dependence of the present quasi-relativistic cal- 
culations, completely different shape (ECPl) [79| and en- 
ergy (ECP2) [801 consistent basis sets available for the 
Cs atom have been used. The orig inal spin-averaged 
Gaussian basis sets from Refs [tI, doj were extended by 
additional diffuse and polarization functions while the 
relevant spin-orbit Gaussian sets were directly borrowed 
from the above references. 

The molecular orbitals (MOs) of Cs2 derived by the 
self-consistent field (SCF) method in the D2h point 
group symmetry were then optimized by the solution of 
the state-averaged complete active space SCF (SA-CAS- 
SCF) problem for the lowest {l-5y'^T,+ g, (1-3) ^^^n,^^ 
and {lY'^Au^g states taken with equal weights jSlJ. The 
dynamical correlation effects were introduced by the in- 
ternally contracted multi-reference configuration interac- 
tion method (MR-CI) [t^. The respective CAS consisted 
of the 7(Tu^g, 4:T^u,g and 2Su,g optimized MOs. MR-CI was 
applied for only two valence electrons keeping the rest 
frozen, i.e. in a full two- valence electron CI scheme while 
the ^-independent core-polarization potentials (CPPs) 
with properly adjusted cutoff radii were employed to take 
into account the remaining core-polarization effects im- 
plicitly. The relevant spin-orbit Gaussian basis set coef- 
ficients were scaled in order to reproduce the experimen- 
tal fine-structure splitting of the lowest excited Cs(6^P) 
state [zH- 

To elucidate the impact of the electron correlation ef- 
fect on the present results, the CPP-CI energies and wave 
functions have been repeatedly evaluated by means of 
the lower C2v group symmetry. In this case the respec- 
tive CAS was restricted by the lAa and IOtt optimized 
MOs. The results obtained in both D2h and C2v repre- 
sentations are found to be almost identical. Furthermore, 
the energies and matrix elements obtained in the frame- 
work of the same CPP-Cl procedure by using ECPl and 
ECP2 basis sets coincided with each other to within a 
few percent. 

The resulting PECs for all electronic states of u- 
symmetry converging to the first, second, and third dis- 
sociation limits are depicted in Fig. [5^, while some of 
the relevant diagonal and off-diagonal spin-orbit coupling 
functions are compared in Fig. [51 



B. Multi- Reference Restricted Active Space 
Configuration Interaction (MR-RAS-CI) 
calculations (Temple University, by author SK) 

At Temple University, we have performed non- 
relativistic as well as a relativistic electronic structure 
calculation to determine the strength of the spin-orbit 
coupling between the ^"'^Sj and b^Ilu potentials of Cs2. 
One electron occupied orbitals are obtained from an 



13 



atomic Hartree-Fock or Dirac-Fock calculation, respec- 
tively. Virtual, highly-excited orbitals are Sturm-type 
functions. The orbitals are labeled Is, 2s, 2p . . . etc in 
analogy with the principal quantum number and orbital 
angular momentum of the hydrogen atom. A configu- 
ration interaction based on molecular determinants se- 
lected by a multi-reference restricted active space (MR- 
RAS-CI) method is used [s^l. Details on the implemen- 
tation of this method in our calculations of heavy di- 
atomic molecules are given in Ref. [83]- Finally, for the 
non-relativistic CI calculation we have evaluated matrix 
elements of the spin-orbit operator 



Hso ^^^^ ^l^-s^ 



N 



iN 



1 I- - 



(s, + 2s-,) , (10) 



where a is the fine structure constant and Si is the spin of 
electron i. The first term of Eq. (ITO| is an one-electron 
operator that describes spin-orbit interactions between 
one nucleus and one electron. Here riN is the separation 
between the i-th electron and nucleus N with charge Zj^ 
and li is the electron orbital angular momentum rela- 
tive to nucleus N. The second term of Eq. ([T0|) is a 
two-electron operator describing the spin-orbit interac- 
tion between electrons. Here fij is the separation be- 
tween electrons i and j and pi is the momentum of elec- 
tron i. 

The closed or filled orbitals up to the Ad shell of Cs 
form the core orbitals of the molecular determinants used 
in the CI for both the non-relativistic and relativistic 
calculation. No excitations are allowed from these shells. 
The 5s^ and shells are core- valence orbitals and are in 
the active space from which single and double excitations 
are allowed. The 6s and 6p orbitals are also added to the 
active space and single, double, and triple occupancy is 
allowed. Finally, we use four each of the s,p,d, and / 
virtual Sturm orbitals to complete the active space. Up 
to double occupancy is allowed for these virtual orbitals. 

The off-diagonal spin-orbit matrix elements of the non- 
relativistic calculation is shown in Fig. [S] We find that 
the one-electron spin-orbit operator of Eq. ^TU\\ provides 
~ 99 % of the total value of the spin-orbit coupling. 

The relativistic potential calculations give an avoided 
crossing between A and b potentials. The smallest en- 
ergy difference AE is evaluated at r =5.5 A and equals 
326 cm^^. According to degenerate perturbtion theory, 
half of this value is equal to the off-diagonal ^^°q spin- 
orbit function at the same internuclear distance. Com- 
parison of ^^°{,Q at r =5.5 A obtained by non-relativistic 
and relativistic calculations shows a good agreement. To 
be precise, Ca6o(5-5 A) = 163.0, 161.6, 158.2 and 164.4 
cm-\ from MR-RAS-CI, ECPl, ECP2 and experiment, 
respectively. 



V. TRANSITION AMPLITUDES TO X STATE 
LEVELS 

The distribution of intensities for transitions from or 
to the X^S^ state is important in testing model poten- 
tials and wavefunctions, for designing data acquisition 
procedures, and is also important for designing transi- 
tion sequences to produce cold Cs2 molecules from cold 
Cs atoms or from Cs Feshbach resonance states. Transi- 
tion intensities can be evaluated by using non-adiabatic 
wavefunctions of the A ^ b complex calculated by FD or 
by DVR methods, and this section will present examples 
of each. 

Radiative properties of the A ^ b complex. To 

test the reliability of the deperturbation analysis de- 
scribed above, we have evaluated the A ~ — >■ X tran- 
sition probabilities, /, and radiative lifetimes r of the 
A ^ b complex according to the relations 



1 



= — y 

TAr^b 2?ic 

VX 



x\{(l>A\dAx\vx)\ 

\{^A\dAx\vx)\^ 



(11) 



A^b^Xl 



where i^A^b^x = E^^-^^iJ') - Ex{vx;Jx) is the 
wavenumber of the rovibronic ~ ^'^11+ — > X^'^'^ 

transition. dAxi^) is the ab initio spin-allowed A^Y,^ — 
X^S^ transition dipole moment calculated in the present 
work by the ECP-CPP-CI method. £;™'=(J') is the 
energy and \4>a) is the non-adiabatic wave- function of 
the A-state component which were obtained by the FD 
method from numerical solution of close coupled radial 
equations with the present empirical PECs and SO func- 
tions. The rovibronic eigenvalues Ex and eigenf unctions 
of the A"^S+ state were obtained by solving the sin- 
gle channel radial equation using the empirical potential 
from [53| . 

The Riga group obtained FTS LIE fluorescence inten- 
sities over a range of X state vibrational levels from a 
J' =238 level with predominantly (calculated to be 88%) 
VA — 'i character. Figure [T51 shows the two predominant 
wavefunction components of this state together with re- 
sults of two computational results for the ab initio tran- 
sition dipole moment. The inset shows this moment over 
an extended range of r. The relative fluorescence intensi- 
ties in Figure [HI show good agreement between observed 
and calculated values. 

The TAr^b — 26 /iS value predicted for the metastable 
ground Ub = level of the fo^H^f^ state, which has only 
0.5% admixture of the singlet A state, is very close to 
the experimental decay rate < 5 x 10** s~^ measured by 
time-resolved fluorescence spectra of Cso dimer molecules 
immersed in a solid helium matrix [84| . The TA^b = 
0.5 /iS obtained for the excited — 37 level, which has 
25% fraction of A^Y,^ character, also agrees well with 
the experimental decay probability of 2.5x 10^ s""'^ of the 



X^E+ transition estimated from the kinetics of 



Cs2 fluorescence measured as a function of temperature 
and Xe buffer gas density in 
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FIG. 13. (Left axis)Wavefunction components for an excited 
state level of J=238, which consists mostly of A''T,^{v = 2), 
but also with a b^Tlg^ component that is shown with dashed 
lines, indicated by two arrows, (right axis) Calculated tran- 
sition dipole moment between the A^'Ei and X^T,'^ state, 
in a.u. The ECPl and ECP2 results from the ECP-CPP-CI 
method are shown as discussed in the text. The inset shows 
theoretical predictions for the transition dipole moment out 
to larger values of r. 



Transition intensities. It was noted that the inten- 
sities for OODR polarization spectroscopy experiments 
at Temple University were much weaker than for anal- 
ogous experiments with Rb2. An explanation for this 
was obtained from Franck-Condon (FC) factors, which 
showed that over the spectral region of interest here, the 
FC factors were about 100 times smaller than the FC fac- 
tors that pertained to the observations of Rb2 A/b ^ X 
transitions reported in [i^l- A plot of FC factors to ^4/6 
levels of interest from X state levels with < 6 is shown 
in Fig. [13 for J— 50. Fluctuations in these FC fac- 
tors are due to the varying fraction of A^T,:^ character. 
Analogous plots for different values of J are qualitatively 
similar once the energy scale is adjusted for the effects of 
centrifugal distortion. 

Note in Fig. (THlthat there are reasonably favorable FC 
factors from low levels of the X^'E'^^ state to low levels 

of fc^IIo^. For the X state, Tg = 4.645 A, while for the 
60 state, Ve — 4.457 AjH^I- Because the values are 
similar, the v=0 wavefunctions overlap. Spin-orbit (CAbo) 



mixes in some A state character into the v=0 of b^IlQ^^, 
producing a significant {bO~ X){0, 0) FC factor, as shown 
in Fig. [151 For higher vibrational levels of b^IlQ^, the 
overlap with X(v — 0) is less. This pattern is repeated for 
higher v of the X state, but extending over more v levels 
of 63n+„. The FC factors increase again when the upper 
state energy approaches v=0 of the state. Thus in 

Fig. [T31 there are in general two regions of appreciable 
FC factors, with a gap between them. 

As noted above, one remarkable aspect of the transi- 



FIG. 14. Relative intensities for the fluorescence of the 
J' = 238 excited state, whose wavefunction components are 
shown in the previous figure, to vibrational levels of the X^E^ 
state, as measured in Riga (wider, open bars) and as calcu- 
lated by the FD method (narrower, filled bars). Intensities 
are normalized to a maximum value of unity in each case. 
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FIG. 15. Logio of Franck-Condon factors for transitions be- 
tween several levels of the X^S^ state to mixed A/b levels, 
for J— 50. Only the A^T,i part of the mixed wavefunction 
was used in the calculation. 



tions observed in Innsbruck from X{v ^ 73) to A/b state 
levels was that the upper level of the strongest transi- 
tions often was primarily triplet in character, as judged 
by the calculated rotational structure or the fractional 
composition. This behavior was predicted by ab ini- 
tio calculations by N. Bouloufa and O. Dulieu, as well 
as by calculations with fitted potentials and DVR wave- 
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FIG. 16. Logio of Franck-Condon factors for transitions be- 
tween two nearby A/b levels to a range of levels of the X^'E'^ 
state, showing contrasting behavior. The inset shows that the 
lower state has a slope with J( J + 1) characteristic of mostly 
b'^Tlou+ character, while the level at 10074.8 cm~^ has a slope 
with J(J + 1) characteristic of mostly A^T,^ character. 
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FIG. 17. Wavefunctions for the excited states in the previous 
figure, together with wavefunctions for the X state levels that 
give the maximum Franck-Condon overlap with the A state 
component. Although the level at 10049.7 cm~^ is mostly 
&"^nou+ character, it is intermixed with A state character to 
give a substantial Franck-Condon factor for a transition from 
X{v = 45). 



functions. Figure [16] shows the contrasting behavior in 
transitions between X state levels and mixed A/b lev- 
els, depending on whether A or 6 character is dominant. 
Figure [T71 shows that a small admixture of char- 
acter, with the correct phase, can produce a significant 
Franck-Condon overlap even for a level that is predom- 
inantly 6^1X011+ in character. Thus the levels that have 
predominantly b state character are preferred as inter- 
mediate levels in the 4-photon STIRAP scheme jT2l - [l3 | 
because they give a much stronger transition on the first 
leg of the second two-photon transition and a more bal- 
anced distribution of transition strengths than levels of 
predominantly A character. 



VI. CAN HYPERFINE STRUCTURE BE 
OBSERVED? 

Hyperfine structure (hfs) was not observed in the ex- 
periments reported here, even though A and b states 
of Na2 exhibited hyperfine structure of several hundred 
MHz (in a regime in which the fraction of fe"^n2u 
character was significant), and the 6 '^S Fermi contact 
splitting in Cs is 9192.63 MHz, as compared with the Na 
3 hfs splitting of 1771.616 MHz Since the OODR 
polarization technique with narrowband lasers is inher- 
ently Doppler-free, one might ask why hfs has not been 
observed in the present work. Furthermore, the capabil- 
ity of exciting alkali dimers from ultracold ground states 
or from Feshbach resonances, either of which are pro- 



duced from cooled atoms, should also lead to possibilities 
for observing hyperfine structure with negligible Dopplcr 
width, in the lowest excited states of these species. With 
these questions in mind, we present a review of available 
information and then make some estimates of possible 
magnitudes of hyperfine structure. 



A. Hamiltonian Matrix Elements 

A detailed theory of hfs in molecular '^H states was 
developed in i8l\ . and applied to molecular beam obser- 
vations of the a^n state of "CO [H,!!!. [H translated 
the parameters of (87| into the more common notation of 
9C I . (soj corrected errors in the Hamiltonian matrix of 
I and also in the forms adapted from [g^l ■ An alterna- 
tive formalism developed by |9l| . originally intended for 
application to the I2 molecule, was applied to the b^Hu 
state of Na2 by Later, [13 



I'^Ag state of Na2, and 



from observations of the 
from data on mixed A^T,^ 



and b^Hu states of Na2, found that their observations 
could be explained almost entirely by the Fermi-contact 
parameter alone. Other studies of hfs in alkali dimers, 
with one exception noted below, have found the Fermi 
contact parameter sufficient. However, we will present 
the Hamiltonian matrix elements for all the magnetic 
dipole parameters, neglecting possible contributions from 
electric quadrupole terms. 

In terms of the hyperfine parameters defined in [9^ . 
(similar to the parameters of [93], but with bp = b + 
c/3), plus two more parameters, d and e, not defined in 
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these references, the elements of the magnetic dipole hfs 
operator are 



{^^Qu±\Hmd\^^1u} 

(^I1iu\Hmd\^^1u) 

{^T1iu\Hmd\^^2u) 
{'^U2u\Hmd\^^2u) 

where 

x^J{J+l); 
F = J + I. 



(2//3)(x/2)i/2(cTd-36F) 
ya 

iy/3)[ix-2)/2]'/'ic-3bF) 
2y{a + bF + 2c/3) (12) 



[i^(F+l)- J(J+l)-/(/+l)], 



2J(J + 1) 



(13) 



a, bp and c are the nuclear spin-orbital angular momen- 
tum interaction, the Fermi-contact interaction, and the 
electron spin-nuclear spin dipolar interaction, d and e are 
two additional dipole interaction terms, proportional, re- 
spectively, to Dii and Dqi in |93i] and [37,]. 

Below, following we define the parameters used 
above and, in the second equality, translate them to 
the parameters used by [37] and [1^. Here C = 
ffS5wMsAiw(Mo/47r), where 5s(mb) and gN{^J'N) are the 
electron and nuclear g-factors (Bohr magneton and nu- 
clear magneton), respectively, and fio is the vacuum per- 
meability. Sums below are over the electrons, i, and nu- 
clei, a. 

» = Y.'-^/ssKIA = i\l/rlJA = 1) = G„/2 



V^E<(|)<a^W-..)|a^i>-(^)- 

^ - E ^C(A = l|(3cos^ 01^ - l)/r!jA = 1) = (^0L) 



^ = E oC(A - -l|e-^'^- sin^ 0,,JrlJA = 1) 



8V15 



(14) 



For e, the expression in terms of 0i^a and r^^Q, while 
straightforward, is too lengthy to give here. 



B. Estimate of Splitting from the Fermi Contact 
and Other Terms 

If we write the Fermi contact term as HpmC = bpl ■ S, 
then from 9l| and [93^, one obtains bp ^ Ahf,atom/4:, 
where Ahf^atom is the atomic ^S* hyperfine parameter. 
For the Cs state, Ahj^atom = 2298.16 MHz (sef. 
Hence the molecular parameter, bp ~ 574.5 MHz. 

As in [33, we will designate a mixed 0+ state by 



where 0, 1, 2 denote '^Ho^„ states with J7 = 0, 1,2, and 
S denotes the A^T,^ state. 

Let us consider a state that is nominally ^Ho„. Consid- 
ering the Hamiltonian matrix elements in Eqs. [21 there 
will be some degree of '^Hi admixture, which we can es- 
timate as 



Cl 



/2xB 



/2xB 



i?(3Hi) - i?(3Ho) 



(16) 



610 



Considering only terms in bp and letting F = F{F + 1) 
etc., the hyperfine shifts of this nominal '^Hq state become 



xyBbf 



[F~J~i]Bbp 



Otso 



fclO 



2Q[F -J- /]kHz, 



(17) 



assuming S^lxlO cm and ^^^q 



' 150 cm-^ For 

J=0, the result is 0. In general, F — |/ — J|...(/ + J). For 
given J, the maximum value of the expression in brackets 
is 2J/, and thus for 1 = 7 (the maximum / value for 
CS2), the maximum value of this term is 280 J kHz. In 
reality, most ^no„+ state levels of interest will be mixed 
with so this estimate should be taken as an upper 

limit. 

Thus, although splitting of the Cs atomic 6 ^S* state 
is 9.2 GHz, hyperfine shifts from the Fermi contact term 
in Cs2 0+ A/b states of low J are less than one MHz, 
due to the relatively large fine structure splitting in Cs2. 
Alternatively, one can say that for 51 = and ft = 1, the 
electron spin precesses rapidly around the internuclear 
axis, thereby reducing terms in I ■ S. 

To our knowledge, only one other hfs term besides the 
Fermi contact term has been convincingly determined in 
alkali dimer states and that is the c term, which was 
found to vary from zero to 9% 6^- in various vibrational 
levels of the c^T,+ state of NaK [H (the many fitted 



|0^) =co|0)-fci|l>+C2|2>+cs|S>, 



(15) 



parameters of 37[ must be reviewed in light of the con- 
clusions of [93*]). From Eqs. (|T2|) . the c term would af- 
fect the observed hfs shifts to a fraction comparable to 
that reported by [9^. If '^Hi„ levels of Cs2 happen to 
be observed in cold molecule spectroscopy, then because 
the ('^Hi„|i?MDpni„) element is diagonal, the term in 
a might be detected. If a were also ^ 5% of bp, shifts 
could be several MHz. 



VII. CONCLUSIONS AND 
ACKNOWLEDGMENTS 

In the interest of providing a model of the A and b 
state potentials and energy level structure for use in the 
experiments directed to the production of cold molecules, 
we have made use of data from several sources, obtained 
for various purposes, as noted above. Because the data 
are relatively sparse in certain regions, the analysis has 
required special procedures. 

Questions that warrant further study have been noted 
above. For example, it appeared that some of the 6"^Hou- 
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(/ parity) low vibrational levels, from the data of [sij . 
could better be assigned to '^Il2u- As laser techniques im- 
prove, it would be interesting to attempt to access these 
levels with high resolution techniques. Secondly, the 
second-order spin-orbit (S02) corrections to the b^Hiu 
potential are quite large, and sensitively dependent on 
the relative potentials of perturbing states, as shown in 
Fig. [5j More extensive data on b^Iliu levels might shed 
light on these S02 corrections. And finally, as discussed 
in Sec. IVI[ it remains a challenge to observe hyperfine 
structure in levels with substantial b^Il2u character. This 
also would require enhanced laser techniques, such as 
used recently to study B^Ilu transitions in Cs2 
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